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ABSTRACT 


This  report  describes  the  modeling  of  measurement  errors  in  precision, 
long  range  tracking  radars  due  to  the  combination  of  environmental  effects 
and  variations  in  target  cross  section.  The  environmental  effects  considered 
are  tropospheric  and  ionospheric  refraction  and  time  delay,  ionospheric 
scintillation,  and  multipath  through  the  sidelobes  at  low  tracking  angles. 

In  addition,  radar  errors  due  to  thermal  and  galactic  noise,  equipment  jitter, 
and  bias  are  included.  Six-degree-of-freedom  vehicle  dynamics  are  modeled 
to  realistically  represent  the  changing  radar  cross  section  due  to  target 
dynamics  such  as  tumbling  and  precession. 
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1 . 0 INTRODUCTION 


ARIES  is  a system  simulation  computer  program  developed  by  Lincoln 
Laboratory  to  study  radar  tracking  and  command  guided  intercepts  in  a 
realistic  radar  environment.  Written  in  FORTRAN  and  designed  for  execution 
on  the  CDC  6600  computer,  it  has  considerable  versatility  in  the  specifica- 
tion of  radar,  target,  tracking  and  environmental  models. 

1 . 1 Program  Purpose 

ARIES  is  designed  to  be  a useful  analytical  tool  for  several  allied 
areas.  Originally  ARIES  was  used  in  the  strategic  BMD  area  to  estimate  the 
metric  state  vector  (position  and  velocity)  of  tracked  targets,  and  then  to 
extrapolate  ahead  in  time  to  determine  an  intercept  point  for  an  ICBM.  The 
radar  measurements  were  subject  to  environmental  effects  which  were  reflected 
in  the  intercept  miss  distances.  Refraction  and  scintillation  models  were 
used  in  ARIES,  and  the  effects  of  various  calibration  schemes  on  target 
location  accuracies  were  studied.  ARIES  was  also  used  to  study  the  problem 
of  multipath  in  low  angle  tracking  and  to  examine  the  effectiveness  of 
various  proposed  schemes  to  overcome  degradations  in  prediction  caused  by 
multipath.  The  use  of  ARIES  in  BMD  studies  was  terminated  in  the  summer  of 
1974. 
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More  recently,  a modified  endoatmospheric  version  of  ARIES,  known  as  the 
HWLPTR  (HOWLS  Projectile  Tracking  Radar)  Program,  has  been  used  in  the  tacti- 
cal area  for  the  evaluation  of  a radar's  performance  in  "backtracking"  an 
incoming  artillery  or  mortar  shell  to  determine  its  point  of  launch.  The 
resulting  point  estimation  error  CEP  values  also  assist  in  the  evaluation  of 
the  drag  model  error  effects  and  in  the  study  of  the  overall  performance  of 
hostile  weapon  location  systems. 

1 . 2 Program  Features 

Since  the  ARIES  Program  provides  a fairly  elaborate  system  simulation, 
it  is  useful  to  tabulate  the  various  features  incorporated  in  ARIES.  The 
major  components  of  the  simulation  are  summarized  below: 

1.  Target  trajectories — accepts  input  of  target  state  vectors  in 
several  coordinate  systems.  Trajectories  with  launch  angle,  reentry 
angle  or  minimum  energy  constraints  from  a given  launch  location  to 
an  impact  point  are  also  available. 

2.  Radar  models — both  mechanically  steered  (dish)  and  phased  array 
radars  are  modeled.  Radar  sensitivity,  beamwidth,  frequency  and 
location  are  specified  by  inputs.  Range  and  angle  measurement 
precision  are  also  specified  by  inputs. 

3.  Radar  measurement  modeling 

a.  Target  modeling — static  cross-section  measurements  on  real 

targets  are  used  in  conjunction  with  rigid  body  dynamics  (Euler's 
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equations)  to  obtain  realistic  dynamic  RCS  simulations.  Constant 
and  sinusoidal,  as  well  as  an  analytic  cylinder,  RCS  options  are 
also  available. 

b.  Noise  and  propagation  effects — radar  measurements  are  corrupted 
by  receiver  noise  (S/N  dependent)  range  independent  noise  effects 
(jitter,  quantization,  etc.)  and  uncorrected  propagation  effects. 
Tropospheric  and  ionospheric  refraction  are  assumed  to  be  cor- 
rected to  within  random  percentages  (input  parameters)  of  the 
true  values.  Ionospheric  scintillation  and  multipath  effects 
corrupt  the  data  but  are  assumed  to  be  uncorrected  in  the 
measurement  model . 

4.  Trajectory  estimation  (Target  Tracking) — Maximum  Likelihood  Estima- 
tion (MLE)  of  the  target  trajectory'  is  performed  based  on  measurement 
data  collected  at  specified  PRF's  over  specified  track  intervals. 
Individual  measurements  are  weighted  according  to  their  measurement 
variances.  ARIES  could  be  easily  extended  to  use  recursive  tracking 
algorithms. 

5.  Target  Discrimination — (Not  presently  implemented.)  Conceptually, 
discrimination  algorithms  would  be  implemented  to  determine  whether 
a particular  simulated  target  constituted  a threat  to  the  defended 
area. 

6.  Interceptor  modeling — (Not  presently  implemented.)  Flight  charac- 
teristics of  one  or  more  interceptor  types  would  be  utilized  to 
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conduct  a command  guided  intercept.  Currently,  the  program  extra- 
polates the  estimated  and  true  target  state  vectors  to  various  time 
(or  altitude)  points  after  termination  of  track  to  obtain  miss 
distances.  Miss  distance  statistics  are  computed  from  the  accumulated 
miss  distances  observed  on  a series  of  Monte  Carlo  tests. 

In  addition  to  the  above  simulation  components,  ARIES  also  accommodates 
multiple  radars,  multiple  targets  and  multiple  track  intervals  on  a given 
simulation  run.  The  feature  of  making  many  Monte  Carlo  runs  for  a given 
scenario  permits  the  generation  of  meaningful  miss  distance  statistics.  A 
building  block/ subroutine  program  structure  lends  itself  to  reasonable 
straightforward  modifications  of  or  additions  to  the  program. 

The  input/output  of  the  ARIES  Program  is  engineer  oriented.  For  input, 
simulation  data  cards  are  conveniently  grouped  into  "packets"  (each  packet 
defines  a target  model,  a radar  model,  an  environmental  model,  etc.)  which 
the  engineer  may  simply  stack  up,  together  with  packets  specifying  the  desired 
simulation  "scenario".  For  output,  an  8y  x 11"  ARIES  Test  Report  (see 
Appendix  II  of  Reference  1)  is  generated  which  provides  the  engineer  with 
descriptions  of  his  input  model  parameters  and  scenarios,  along  with  the 
resultant  simulation  data  and  statistics.  The  outputs  are  all  organized 
into  logical  sections  which  are  indexed  for  ready  reference.  Outputs  from 
ARIES  also  include  trajectory  plots  superimposed  on  a world  map,  plots  of 
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true  and  measured  target  cross-section,  and  a radar  measurements  tape  con- 
taining metric  and  RCS  data  for  processing  by  other  programs. 

1 . 3 Program  Documentation 

The  ARIES  Program  is  documented  in  three  separate  Lincoln  Laboratory 
Technical  Notes  as  follows: 

1 . The  ARIES  Program  - A General  Overview  and  Users ' Guide 

2.  The  ARIES  Program  - Coordinates,  Transformations,  Trajectories  and 
Tracking 

3.  The  ARIES  Program  - Analysis  and  Generation  of  Simulated  Radar 
Measurements 

The  first  report  presents  a general  discussion  of  the  ARIES  Program, 
including  the  logical  organization  of  the  program  and  descriptions  of  all 
subroutines.  All  of  the  options  available  to  a user  are  discussed  and  the 
methods  of  setting  up  the  input  "packets",  including  controls  to  activate  the 
various  options,  are  presented.  A typical  run,  including  a complete  ARIES 
Test  Report  (output),  is  discussed  in  this  first  volume. 

The  second  and  third  volumes  contain  all  of  the  relevant  mathematics 
and  the  models  used  in  TRIES.  Most  of  the  deterministic  mathematics  (co- 
ordinate systems  and  tr;msformations , trajectory  generation  and  estimation, 
miss  distance  calculations,  etc.)  are  in  the  second  report.  This,  the  third 
report,  is  primarily  concerned  with  the  generation  of  radar  measurements, 
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including  the  corruptive  effects  of  noise,  radar  biases,  propagation  and 
time-varying  radar  cross-section. 

1.4  Organization  of  This  Report 

In  Section  2 we  start  by  examining  the  parameters  and  equations  which 
are  used  to  generate  "measured"  radar  data  for  a Monte  Carlo  analysis. 
Equations  are  presented  for  both  a dish  type  and  phased  array  radar.  Section 
3 deals  in  more  detail  with  random  noise  errors  (thermal,  galactic,  etc.)  and 
bias  errors  (quantization,  etc.),  and  includes  a description  of  the  parameters 
needed  to  determine  the  signal  to  noise  ratio  (SNR) . In  Section  4 a model 
is  presented  for  those  types  of  atmospheric  refraction  which  most  severely 
affect  radar  performance.  Both  ionospheric  and  tropospheric  refractions  are 
considered  and  the  proper  model  characteristics  are  provided  for  both  day 
and  nighttime  use.  It  is  also  assumed  that  a given  percentage  of  the  refrac- 
tion error  may  be  corrected  and  a parameter  to  specify  the  amount  of  this 
residual  error  is  provided.  Section  5 is  concerned  with  the  effects  of 
ionospheric  scintillation  on  the  signal  to  noise  ratio  and  radar  angle  errors. 
A model  is  presented  based  on  measurements  made  at  M.I.T.  Lincoln  Laboratory 
Millstone  Radar  and  earlier  work  by  Nakagami  [13].  In  Section  6 multipath 
errors  due  to  a second  signal  path  reflected  into  the  sidelohes  are  examined. 
Typical  antenna  sidelobe  patterns  and  reflection  coefficients  are  discussed 
along  with  methods  for  correcting  the  model  for  different  mainlobe  beamwidths. 
The  last  section,  Section  7,  is  concerned  with  a description  of  the  RCS 
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computations  made  by  the  program.  Several  RCS  calculation  techniques  are 
presented  including  the  use  of  stored  static  patterns,  and  purely  mathematical 
models  of  some  typical  object  shapes. 
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2 . 0 RADAR  MEASUREMENTS 


The  ARIES  Program  makes  use  of  a Monte  Carlo  technique  to  analyze  the 
errors  produced  by  a radar  system.  However,  the  use  of  a Monte  Carlo  method 
requires  the  ability  to  accurately  simulate  the  radar  system  under  test.  In 
particular,  it  must  be  possible  to  simulate  radar  measurement  errors  at  any 
point  along  a target  trajectory.  These  measurement  errors  are  obtained  by 
corrupting  the  true  coordinates  of  the  target  by  factors  depending  upon  the 
radar  system  parameters,  signal-to-noise  ratio  and  propagation  effects.  If 
the  radar  specified  is  an  RAE  tracker  (dish  type)  then  these  measurements  are 
supplied  in  RAE  coordinates.  For  a phased  array  or  RUV  tracker,  then  the 
measurements  are  returned  in  RUV  coordinates.  In  the  following  two  sections 
the  equations  used  to  corrupt  the  true  radar  coordinates  will  be  presented  for 
the  RAE  and  RUV  cases.  Detailed  explanations  of  each  term  in  the  measurement 
equations  may  be  found  in  the  later  sections  which  describe  the  development  of 
the  noise  and  environmental  models.  We  start  with  the  RAE  tracker. 

2 . 1 Radar  Measurements  in  RAE  Coordinates 
th 

For  the  i—  trajectory  point  at  which  a measurement  is  desired,  let  Rj,(i) , 
Ap(i)  and  E^(i)  represent  the  true  range,  azimuth  and  elevation.  Then  denote 
by  R^(i,j),  A^(i,j)  and  E^(i,j)  the  measured  values  of  range,  azimuth  and  ele- 
vation of  the  i—  trajectory  point  on  the  j—  Monte  Carlo  run.  Note  that 
although  the  true  values  for  the  RAE  coordinates  are  the  same  on  each  Monte 
Carlo  run,  the  measured  values  will  not  be,  hence  the  need  for  a variable  j to 
index  the  Monte  Carlo  run  number. 
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The  equations  for  the  R,  A,  and  E coordinate  measurements  are: 


fyq  ( i > 3 ) ~ ( i ) 

+ kl  °RJ 

+ k2  aRT^^  ^2^’^ 
+ k3  ABR^-^ 

+ k4  oRI(i)  RN3C3) 

+ k5  aRTp,1-)  ^4*-^ 

= Aj.fi, j) 

+ k.j  a^j  RNg(i»j) 

+ k2  oAT(i)  RN6(i,j) 

+ k3  aBa^) 

+ k^  AASint 

Ejq  ( i > j ) = ( i > j ) 


Jitter 

Thermal 

Bias 

Ionospheric  refraction 
Tropospheric  refraction 

Jitter 

Thermal 

Bias 

Ionospheric  scintillation 

Jitter 

Thermal 

Bias 

Ionospheric  refraction 
Tropospheric  refraction 
Multipath 


+ kl  °EJ  ^n<7  C i > j ) 

+ k2  a£T(i)  RNg(i,j) 

+ k3  ABE^ 

+ ^4  aEI  ^ ^Ng(j) 

+ k5  aETp(-1')  ™10W 
+ k6  AEMul 

where : 

T denotes  the  true  or  uncorrupted  coordinates 


(2-1) 


(2-2) 


(2-3) 
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M denotes  the  measured  or  simulated  coordinates 
i is  an  index  for  the  track  points 
j is  an  index  for  the  Monte  Carlo  run  number 

k-^-ky  are  switches  which  are  1 if  the  particular  form  of  corruption 
(thermal,  bias,  etc.)  is  to  be  used,  0 if  it  is  to  be  bypassed. 

RN  (i,j)  are  independent  samples  from  a zero  mean  unit  variance  Gaussian 
random  number  generator.  The  notation  RNn(i,j)  is  used  to  indi- 
cate that  a new  sample  from  the  random  number  generator  is 
required  for  each  track  point  of  each  Monte  Carlo  run. 

RNn(j)  are  independent  samples  frcm  a zero  mean  unit  variance  Gaussian 

random  number  generator.  The  notation  RNn(j)  is  used  to  indicate 
that  the  random  number  sample  remains  constant  over  the  track 
points  and  changes  only  at  the  beginning  of  each  Monte  Carlo  run. 

The  remaining  parameters  are  described  in  sections  3.0  through  6.0.  The 
table  below  contains  a list  of  these  parameters  along  with  the  section  number 
where  a full  description  is  available. 


Parameter 

Descr iption 

Section  ] 

°RJ  * °AJ’  °EJ 

Jitter  errors 

3.1.1 

°RT  * aAT’  °ET 

Thermal  errors 

3.2.1 

ABr  , AB^ , AB^ 

Bias  errors 

3.3.1 

VALn(j) 

Constant  or  random  sample  used  with 
bias  errors 

3.3.1 

°RTp ’ °ETp 

Residual  tropospheric  refraction  errors 

4.3.1 

°RI ’ °EI 

Residual  ionospheric  refraction  errors 

4.3.2 

AASint 

Ionospheric  scintillation  angle  error 

5.4 

AEMul 

Multipath  elevation  error 

6.2 
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2.2  Radar  Measurements  in  RUV  Coordinates 


For  a phased  array  radar  the  measurement  simulation  is  performed  in  RUY 
coordinates.  Again,  let  i index  the  trajectory  points  at  which  a measurement 
is  desired,  and  Rju) , lij(i) , \j(i)  denote  the  true  values  of  the  RUY  coordin- 
ates at  point  i.  The  measured  values  are  represented  by  R. j( i , j ) , U^(i,j), 
Vv(i,i)  where  j is,  as  in  the  last  section,  the  Monte  Carlo  run  number. 

The  measurement  equations  for  the  R,  U and  V coordinates  are  of  a similar 
fora  to  equations  (2-1)  to  (2-3)  and  are  as  follows: 


RvjliJ)  ■ Rj-(i) 

+ kl  :RJ 

+ k2  '"RT1-1' 

+ k5  AB^j) 

+ k4  -"ri^  ^3(2 ) 

+ k5  rRTp(l1  ^4  ^ 


(2-4) 

Jitter 

Thermal 

Bias 

Ionospheric  refraction 
Tropospheric  refraction 


UM(i,i)  = UT(i,i) 

+ k1  -yj  R-N'5(i,j) 

+ k-,  ~LT(i)  RX6(i,j) 
- k5  ABydJ) 

+ k4  rUTp 
+ k5  Jyj  RN3O) 

+ k6  AUMul 
+ k7  -USint 


C-5) 

Jitter 

Thermal 

Bias 

Tropospheric  refraction 
Ionospheric  refraction 
Multipath 

Ionospheric  scintillation 
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VM(i,j)  = VT(i,j) 


(2-6) 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


kl  °VJ  RN9^i»^ 

k2  °vt^  mioV>n 

k3  ABvrCi ,j) 
k4  °VTp  ^11^ 
k5  °VI  ^12  ^ 
k6  AVMul 
k7  AVSint 


Jitter 

Thermal 

Bias 

Tropospheric  refraction 
Ionospheric  refraction 
Multipath 

Ionospheric  scintillation 


where : 

T denotes  the  true  or  uncorrupted  coordinates 
M denotes  the  measured  or  simulated  coordinates 
i is  an  index  for  the  track  points 
j is  an  index  for  the  Monte  Carlo  run  number 

kj -ky  are  switches  which  are  1 if  the  particular  form  of  corruption 
(thermal,  bias,  etc.)  is  to  be  used,  0 if  it  is  to  be  bypassed. 

RN  (i,j)  are  independent  samples  from  a zero  mean  unit  variance  Gaussian 
random  number  generator.  The  notation  RNn(i,j)  is  used  to 
indicate  that  a new  sample  from  the  random  number  generator  is 
required  for  each  track  point  of  each  Monte  Carlo  run. 

RN  (j)  are  independent  samples  from  a zero  mean  unit  variance  Gaussian 

random  number  generator.  The  notation  RN’n(j)  is  used  to  indicate 
that  the  random  number  sample  remains  constant  over  the  track 
points  and  changes  only  at  the  beginning  of  each  Monte  Carlo  run. 
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The  remaining  parameters  are  described  in  sections  3.0  through  6.0.  The 
table  below  contains  a list  of  these  parameters  along  with  the  section  number 
where  a full  description  is  available. 

Section  No. 

3.1.2 

3.2.2 

3.3.2 

3.3.2 

4.3.1 

4.3.2 
5.4 

6.2 

2.3  Description  of  Measurement  Generating  Subroutines 

The  measurement  equations  of  the  preceding  section  are  implemented 
predominately  in  subroutine  MEASUR  which  is  flowcharted  in  Figure  2-1.  In 
actual  use,  MEASUR  is  called  by  the  monitor  subroutine  MQNITR  whenever  a 
simulated  radar  measurement  is  desired.  The  unerrored  or  true  coordinates  are 
obtained  and  "corrupted"  by  MEASUR,  and  these  errored  coordinates  are  returned 
to  the  track  file.  Also  the  jitter  one  sigma  value  a^j,  and  thermal  one  sigma 
value  a^(i)  are  combined  to  form  a variance 


Parameter 


7RJ’  aUJ’  aVJ 


aRT’  aUT’  aVT 

ABr,  ABy,  ABy 

VAL  (j) 
nu^ 


aRTp ’ aUTp’  aVTp 

aRI ’ aUI’  aVI 

AUC • . , AVC . 
Sint’  Sint 


AUMul’  AVMu1 


Description 
Jitter  errors 
Thermal  errors 
Bias  errors 

Constant  or  random  sample  used  with  bias 
errors 

Residual  tropospheric  refraction  errors 

Residual  ionospheric  refraction  errors 

Ionospheric  scintillation  angle  errors  in 
U and  V 

Multipath  elevation  errors  in  U and  V 
coordinates 
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18-3  16856; 


FORM  ERRORED  COORDINATES 
BY  CORRUPTING  TRUE  COORDINATES 
(EITHER  RAE  OR  RUV  COORDINATES) 


CALCULATE  ERROR  VARIANCES 
FOR  USE  BY  THE  MLE 


CHECK  FOR  ERRORS 
IN  CALLING  MEASUR 
AND  SELECT  ENVIRONMENTAL 
TABLES 


OBTAIN  TRAJECTORY 
AND  RCS  DATA  FROM 
STORAGE  ARRAY 


CALCULATE  TRAJECTORY 
AND  RCS  DATA 

PLACE  TRAJE 
RCS  DATA  |H 
ARRAY 

CTORY  AND 
ITO  STORAGE 

Fig.  2.1.  Flow  Chart  of  Subroutine  MEASUR. 
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2 

CRest 


(i) 


+ ctrt^ 


which  is  supplied  to  the  maximum  likelihood  estimator  (NILE)  at  each  track 
point  i.  Similar  combined  variances  are  calculated  for  the  A and  E (or  U and 
V)  coordinates  and  passed  to  the  MLE. 

An  option  is  also  provided  to  write  a tape  containing  radar  cross  section 
data  and  is  described  in  more  detail  in  Section  7.2. 

In  order  to  conserve  execution  time  several  techniques  are  used.  First, 
all  multidimensional  arrays  are  equivalenced  to  dummy  one -dimensional  arrays 
so  that  subscript  calculation  time  would  be  minimized.  Second,  an  assembly 
language  Gaussian  random  generator  is  used  instead  of  a Fortran  version  at  a 
savings  in  time  of  1/3  to  1/2.  And  finally,  a storage  mode  option  is  provided 
to  avoid  the  recalculation  of  factors  such  as  the  true  RAF  coordinates,  radar 
cross  section  (RCS) , etc.,  on  each  subsequent  Monte  Carlo  run.  This  mode  can 
be  used  only  in  those  cases  where  the  stored  coordinates  (R,  A,  E,  U,  V,  T,  E, 
RCS,  and  ALT  (altitude))  are  not  varying  from  one  Monte  Carlo  run  to  the  next. 
Also  the  measurements  must  be  required  at  the  same  point  on  the  trajectory 
during  each  Monte  Carlo  run.  Under  these  conditions  a savings  of  approximately 
25?o  is  possible  in  the  overall  execution  time. 

With  the  exception  of  ionospheric  scintillation  and  radar  cross  section 
all  the  errors  are  calculated  within  subroutine  MEASUR.  Ionospheric  scintilla- 
tion errors  are  produced  by  a call  to  subroutine  ISCINT  and  radar  cross  section 
calculations  by  a call  to  subroutine  RCS. 
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3 . 0 RADAR  NOISE  AND  BIAS  PARAMETERS 


The  types  of  errors  modeled  fall  into  two  categories;  random  and  bias. 

The  random  errors  are  due  to  factors  such  as  galactic  and  thermal  noise, 
receiver  noise,  analog  to  digital  converter  jitter,  etc.  These  noises  are 
uncorrelated  from  one  measurement  to  the  next  and  may  or  may  not  have  time 
varying  statistical  properties.  On  the  other  hand,  bias  errors  are  considered 
to  be  those  errors  which  remain  in  the  system  after  calibration  and  are  fixed 
from  measurement  to  measurement . These  include  factors  such  as  inaccurate 
calibration  procedures,  survey  errors,  and  mechanical  problems  such  as  feed- 
horn  sag  in  a dish  array. 

It  is  assumed  that  the  various  sources  of  bias  and  noise  are  statistically 
uncorrelated.  This  is  a reasonable  assumption  for  types  of  errors  described 
in  the  next  sections. 

3.1  Jitter  Errors 

These  are  errors  whose  statistics  do  not  vary  as  a function  of  the  target 
range,  although  the  samples  at  adjacent  track  points  are  considered  to  be 
uncorrelated.  Typical  errors  of  this  type  are  the  quantization  errors  in 
analog  to  digital  converters  following  phase  detectors  or  logarithmic  ampli- 
fiers. Also  quantization  errors  in  devices  such  as  shaft  encoders  for  dish 
antenna  angle  measurements  and  range  counters  fall  in  this  category. 

Because  it  is  not  possible  to  accurately  model  the  statistical  distribu- 
tion of  these  errors,  a Gaussian  choice  is  made  as  a reasonable  approximation. 
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At  each  track  point  the  jitter  errors  are  considered  to  be  zero  mean  with 


2 2 2 2 2 
variances  a^j,  a^j,  a^j  for  RAE  coordinates,  and  variances  a^j,  c^j. 


aV J for 


RUV  coordinates. 


3.1.1  Dish  Radar 

If  a dish  radar  (RAE)  is  selected,  then  o^j,  a^j,  and  a^j  are  input 
parameters  read  by  ARIES.  These  parameters  remain  constant  for  all  track 
points  in  the  simulation.  However,  at  each  track  point  new  random  number 
samples  are  obtained  from  the  distribution  to  calculate  the  jitter  errors. 

3.1.2  Phased  Array  Radar 

For  a phased  array  radar  (RUV) , a^j,  and  Oyj  are  the  jitter  one 
sigma  values  read  in  by  ARIES  during  execution.  These  parameters  remain  con- 
stant for  all  track  points  in  the  simulation.  However,  at  each  track  point 
new  random  number  samples  are  obtained  from  the  distribution  to  calculate  the 
j itter  errors . 

3 . 2  Thermal  Errors 

Thermal  errors  are  used  to  denote  noise  errors  whose  statistics  vary  as  a 
function  of  the  target  range,  and  which  produce  errors  uncorrelated  from  one 
track  point  to  the  next.  Examples  of  this  type  of  errors  are  galactic  noise, 
internal  receiver  noise,  and  clutter.  To  a very  good  approximation  the 
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statistics  of  this  type  of  noise  may  be  represented  as  Gaussian  with  a zero 
mean  and  a given  variance  depending  on  the  source  of  the  noise. 

3.2.1  Dish  Radar 

For  the  dish  type  radar  the  one  sigma  values  for  the  R,  A,  and  E coordi- 
nates are 


OrtU)  = oRT 
OatU)  = °At 
aET^1)  = aET 


1 

\ZSNR1 


1 


1 

v/SNR1 


(3-1) 

(3-2) 

(3-3) 


where  SNR  is  the  signal  to  noise  ratio  defined  in  Section  3.2.3. 
Opp  are  the  reference  thermal  one  sigma  values  defined  at  a range 
below)  and  are  read  by  ARIES  as  input  parameters.  aR^(i)  , 
denote  the  thermal  one  sigma  values  at  track  point  i. 


°£T » °at 

R (see  SNR 
o v 

and  OgT(i) 


3.2.2  Phased  Array  Radar 

For  the  phased  array  radar  the  corresponding  set  of  equations  are 

°RTti) ' 0rt  7k 
0utU)  = 
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1 


where  SNR  is  defined  below.  The  remaining  terms  have  a significance  similar 
to  their  RAE  counterparts  above. 

3.2.3  Signal  to  Noise  Ratio 

The  signal  to  noise  ratio  (SNR)  is  defined  as  the  ratio  of  signal  power 
to  noise  power.  In  ARIES  the  SNR  is  used  to  multiply  the  reference  thermal 
one  sigma  value  by  a factor  which  depends  on  the  target  range,  target  cross 
section,  and  ionospheric  scintillation.  In  addition,  for  a phased  array  radar, 
antenna  gain  reduction  due  to  scanning  the  beam  is  included. 

The  RAE  tracker  has  the  following  expression  for  SNR: 


(3-4) 


where 


R is  the  range  of  the  target 

RCS  is  the  left  channel  radar  cross  section  of  the  target. 


G is  the  attenuation  (or  gain)  in  signal  amplitude  due  to  ionospheric 
scintillation,  (see  Section  5) 


2 

Rq  is  the  reference  range  at  which  the  SNR  equals  one  (0  dB)  for  a 1 
radar  cross  section  and  G=1 . 
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For  the  RUV  tracker  a similar  equation  is  obtained  as: 

SNR  = (j^.)  RCS  G4  cos20  (3-5) 

The  R,  RCS,  and  G terms  have  the  same  interpretation  as  above. 

0 is  the  scan  angle  off  the  phased  array  beam  (0=0  corresponds  to  a 

2 

broadside  beam) . The  cos  0 represents  the  effect  of  antenna  gain 
reduction  due  to  scanning  the  beam. 

Rq  is  the  reference  range  for  a broadside  beam  (0=0)  at  which  the  SNR 
equals  one  (0  dB)  for  a 1 m^  radar  cross  section  and  G=l. 

3 . 3 Bias  Errors 

ARIES  allows  bias  errors  to  be  specified  in  either  of  two  possible 
manners.  A fixed  value  for  the  bias  may  be  introduced  which  remains  constant 
for  all  track  points  in  every  Fbnte  Carlo  run.  Or  a random  number  sample 
may  be  generated  at  the  start  of  each  Monte  Carlo  run,  and  used  as  the  bias 
for  all  track  points  on  that  particular  Monte  Carlo  run.  On  the  next  Monte 
Carlo  run  a new  bias  selection  is  made  from  the  random  number  generator. 

The  first  method  for  specifying  the  bias  is  used  if  it  is  assumed  that 
the  bias  is  constant  for  all  conditions.  Note  that  if  this  is  in  fact  the 
case,  and  this  constant  is  known,  then  it  may  be  subtracted  from  the  received 
data  and  cause  no  measurement  errors. 
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The  second  type  of  bias  errors  results  from  the  ability  to  partially 
correct  the  data  by  a calibration  process.  This  type  of  bias  error  is  assumed 
to  lie  in  a range  of  values  centered  about  the  mean  or  true  bias.  From  a set 
of  measurement  runs  the  mean  can  be  estimated  and  then  subtracted  from  the 
data.  Thus  the  remaining  bias  may  be  modeled  as  a constant  value  sample 
function  from  a zero  mean  process  which  is  added  to  the  measurements  at  each 
track  point.  The  constant  value  for  the  sample  function  is  selected  from  a 
Gaussian  distribution,  since  the  exact  form  of  the  distribution  cannot  be 
determined. 


3.3.1  Dish  Radar 

The  bias  terms  for  a dish  radar  in  R,  A and  E coordinates  are 


t> 

< i. 

" aRB  VAL1^ 

(3-6) 

ABA(j) 

= VAL2(j) 

(3-7) 

ABg(j) 

aEB  VAL3^ 

(3-8) 

where  o^,  a and  are  input  parameters  read  by  ARIES  and 


VALn(j) 


K - for  a constant  bias  term 


RN  (j)  - for  a bias  term  constant  only  for  a given 
Monte  Carlo  run 


where  FlNn(j ) is  an  independent  sample  from  a zero  mean,  unit  variance  Gaussian 
random  number  generator.  The  random  number  sample  is  indexed  by  the  Monte 
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Carlo  run  number,  j,  to  indicate  that  a new  sample  is  desired  each  time  j 
changes . 


3.3.2  Phased  Array  Radar 

For  a phased  array  radar,  the  bias  equation  for  range  is  the  same  as  that 
for  a dish  radar.  The  only  type  of  bias  considered  in  U and  V is  that  due  to 
misalignment  of  the  phased  array  boresight  direction.  The  equations  are 

ABr(i)  = 0rb  VALjCj)  (3-9) 

ABy(i,j)  = a^[cosC^'A)cosE/cosEL]VAL.,(j)  (3-10) 

AByCi.j)  = [am U tanEL+c^W]VAL3(j)  (3-11) 

where  a^g,  a^g,  and  a™  are  the  bias  values  read  by  ARIES  as  input.  VALn(j) 

is  the  same  as  described  above  for  the  dish  radar.  The  bias  errors  at  bore- 
sight  are  a^g  and  a^g.  A and  E are  the  target  azimuth  and  elevation  angles 
and  AZ  and  EL  are  the  azimuth  and  elevation  of  the  phased  array  boresight 
direction.  Also,  W = v/l-U^-V^. 

Note  that  for  the  cases  where  the  random  number  bias  technique  is  used, 
three  separate  and  independent  random  number  samples  are  generated,  one  for 
each  coordinate. 
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4.0  ATMOSPHERIC  REFRACTION 


An  important  consideration  which  must  be  taken  into  account  in  the  formu- 
lation and  design  of  any  radar  system  is  the  degradation  in  performance  re- 
sulting from  propagation  effects  in  the  atmosphere.  In  general,  the  regions 
of  the  atmosphere  which  affect  the  propagation  of  electromagnetic  waves  are 
the  troposphere  and  the  ionosphere.  The  effects  which  are  prevalent  when  radio 
waves  traverse  the  atmosphere  manifest  themselves  in  a variety  of  ways.  In 
this  section,  two  of  the  phenomena  in  the  troposphere  and  ionosphere  which 
most  severely  effect  radar  systems  are  described  as  they  apply  to  the  models 
used  in  the  ARIES  program.  In  later  sections  other  atmospheric  effects  will 
be  discussed. 

4 . 1 Atmospheric  Refractivitv  Model 

The  approach  taken  here  is  summarized  as  follows.  A standard  atmospheric 
refractive  structure  for  the  troposphere  and  ionosphere  is  selected  and  this 
structure  is  assumed  to  be  an  average  of  prevailing  conditions  around  the  radar 
station.  Using  this  atmosphere  and  a ray  tracing  computer  program,  the  mean 
refraction  error  as  a function  of  range  and  elevation  angle  is  determined. 

The  mean  refraction  error  represents  the  contribution  to  the  radar  target  po- 
sition error  produced  by  the  atmosphere  before  any  refraction  correction 
technique  is  applied. 

An  assumption  is  made,  however,  that  some  statistical  method  is  used  to 
correct  for  the  effects  of  refraction.  This  correction  method  leaves  a resid- 
ual error  in  the  radar  position  measurement.  The  resulting  residual  error  is 
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simulated  in  the  ARIES  program  by  generating  a random  variable  whose  standard 
deviation  is  some  fraction  of  the  mean  refraction  error.  In  the  following 
sections,  details  of  the  atmospheric  refractivity  models  are  described  and  the 
method  used  for  generating  the  residual  error  is  presented. 

4.1.1  Tropospheric  Model 

The  Central  Radio  Propagation  Laboratory  (CRPL)  Reference  Atmosphere  - 
1958 [7]  is  the  model  used  by  the  ARIES  Program  for  the  troposphere.  In  this 
model,  the  troposphere  is  stratified  in  three  layers  with  the  refractivity  in 
each  layer  defined  by  the  Eqns.  (4-1),  (4-2),  and  (4-4). 

For  height,  h,  between  h^  s h < h^  + 1 km  the  refractivity  is 


N(h)  = Ng  + (h  - hs)AN 


(4-1) 


where 


AN  = -7.32  exp(5 . 577  x 10'3  N ) 
h^  = surface  height 

N = surface  refractivity  (assumed  320  N-unit  for  this  model) . 


For  h + 1 s h s 9 bn 


N(h)  = Nx  exp[ -c(h  - hs  - 1)  ] 


(4-2) 


where 


N,  = value  of  N at  1 km  above  the  surface 


_ In  Nl/105 
c - —3 r 


" 8 - h ’ 

s 


(4-3) 
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and  for  h > 9 km 


N(h)  = 105  exp [-0. 1424 (h-9) ] . 


(4-4) 


The  above  model  primarily  depends  on  specifying  Ns,  the  radio  refractivity 
at  the  surface  of  the  earth.  The  radio  refractivity  at  the  surface  is  a 
function  of  meteorological  conditions  at  the  surface  and  can  be  represented  by 


where  T is  the  air  temperature  in  degrees  Kelvin;  p,  the  total  pressure  in 
millibars;  and  e,  the  partial  pressure  of  water  vapor  in  millibars.  The  first 
term,  ap/T,  applies  both  to  optical  and  radio  frequencies,  and  is  referred  to 


relationship  required  only  at  radio  frequencies. 

According  to  Kerr  [4] , the  constants  a and  b are  approximately  equal  to 
79°K/mb  and  4800°K,  respectively.  Campen  and  Cole  [5]  quote  74.4°K/mb  and 
4973°K.  Through  consideration  of  various  propagation  experiments,  Smith  and 
Weintraub  [6]  arrive  at  the  values  of  77.6°K/mb  and  4810°K.  The  latter  values 
for  the  constants  a and  b are  the  ones  used  in  the  ARIES  Program. 

The  expression  in  Eqn.  (4-5)  for  the  radio  refractivity  is  independent  of 
frequency  in  the  100-MHz  to  10-GHz  range.  According  to  Smith  and  Weintraub 
[t>] , this  expression  is  considered  to  be  valid  to  0.5%  in  N units  for  frequen- 
cies up  to  30  GHz  and  normally  encountered  ranges  of  temperature,  pressure, 
and  humidity. 


(4-5) 


2 

as  the  dry  term,  while  the  second  term,  abe/T^,  is  the  explicit  water  vapor 
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4.1.2  Ionospheric  Model 


The  refractive  index  of  a medium  containing  free  electrons  is  given  by 
Eccles  and  Larmor  [8]  as 


n 


1 


Ne  e2 


e m (i) 

o 


2 


(4-6) 


where 

3 

- electron  density  in  electrons/m 

-19 

e - electron  charge  1.6  x 10  C 

-31 

m - electron  mass  9.1  x 10  kg 

eQ  - electric  permittivity  of  free  space  8.841941  x 10 
a)  - angular  frequency  of  incident  wave 


This  relationship  is  valid  when  the  earth's  magnetic  field  and  the  collision 
of  electrons  with  ions  and  neutral  particles  are  neglected. 

The  distribution  of  electron  density  with  height  above  the  earth's 
surface,  under  equilibrium  conditions,  is  assumed  to  follow  the  formulation 
by  Chapman  [8]  which  is 


Ne  = Nm  exP  " z " exP(~z))} 


(4-7) 


where  Nm  is  the  electron  density  at  the  maximum  level  of  ionization  and  z is 
the  normalized  height  given  by 


z 


h - *m 
H 


(4-8) 
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The  parameters  for  the  normalized  height  are  h^,  the  height  of  maximum 
ionization  density  and  H,  the  height  of  the  homogeneous  atmosphere  at  a given 
temperature . 

Using  Eqns.  (4-6)  , (4-7)  and  (4-8)  , the  ionospheric  refractive  index  is 
modeled  by  assuming  a daytime  ionosphere  which  consists  of  three  layers  (E,  Fp 
F-,)  and  a nighttime  ionosphere  which  consists  of  two  layers  (E,  F)  . The  values 
for  the  parameters  N , h and  H used  in  the  respective  daytime  and  nighttime 
models  are  given  in  table  4-1. 


Time 

Layer 

H - (km) 

g 

l 

Nj^  - (electrons/m^) 

Daytime 

E 

10 

100 

1.5  x 1011 

F1 

40 

200 

3.0  x 1011 

F2 

50 

300 

12.5  x 1011 

Nighttime 

E 

10 

120 

0.8  x 101U 

F 

45 

250 

4.0  x 1011 

TABLE  4-1 

PARAMETERS  OF  IONOSPHERIC  MODEL 
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4.2  Refraction  Errors 


Using  the  atmospheric  refractivity  model  described  in  Section  4.1,  a ray 
tracing  computer  program,  developed  by  the  General  Electric  Co.,  was  used  to 
determine  radar  range  and  elevation  errors.  The  ray  tracing  program  computes 
the  refractive  bending  which  occurs  as  the  electromagnetic  wave  encounters 
variations  in  atmospheric  refractive  index  along  its  propagation  path.  From 
this  result,  the  radar  range  and  elevation  errors  are  computed  as  a function 
of  elevation  angle. 

Tables  (4-2A),  (4-3A)  and  (4-4A)  give  the  range  error,  er,  and  elevation 
error,  eq,  as  a function  of  elevation  angle  used  in  the  ARIES  Program  for  the 
troposphere,  daytime  ionosphere  and  nighttime  ionosphere,  respectively.  These 
tables,  however,  give  the  measurement  error  assuming  the  electromagnetic  wave 
has  passed  through  all  of  the  refracting  medium.  For  waves  which  pass  through 
only  part  of  the  refracting  medium,  a fractional  adjustment  must  be  applied  in 
order  to  obtain  the  error.  The  error  for  range,  AR,  and  elevation,  AE  is 
given  by 


AR  = er(E)  f1G1) 
AE  = ee(E)  f2(h) 


(4-9) 


where  f^  and  f-,  are  the  fractional  adjustment  terms  given  in  tables  (4-2B) , 
(4-3B)  and  (4-4B) , E is  elevation  angle  and  h is  the  altitude  of  the  target. 
Graphs  of  tables  (4-2),  (4-3)  and  (4-4)  are  given  in  Figs.  4-1  to  4-3. 
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Elevation  (mr) 

er(Elev.)  (m) 

ee(Elev.)  (mr) 

0.000 

107.869 

12.280 

17.453 

69.159 

8.580 

34.907 

49.256 

6.250 

69.813 

30.328 

3.930 

104.720 

21.549 

2.810 

139.626 

16.916 

2.200 

174.533 

13.503 

1.760 

261.799 

9.205 

1.180 

349.066 

6.980 

.870 

523.599 

4.816 

.550 

698.132 

3.962 

.400 

872.665 

3.139 

.270 

1047.198 

2.743 

.180 

1221.730 

2.530 

.100 

1396.263 

2.438 

.050 

1483.530 

2.408 

.030 

TABLE  (4-2A) 
TROPOSPHERIC  ERRORS 
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Altitude  (km) 

(Alt.) 

f2(Alt.) 

0.0 

0.000 

0.000 

1.0 

.140 

.080 

2.0 

.260 

.140 

3.0 

.350 

.200 

6.0 

.590 

.330 

10.0 

.790 

.440 

20.0 

.950 

.600 

30.0 

.970 

.720 

60.0 

.980 

.830 

100.0 

.981 

.890 

300.0 

.990 

.940 

600.0 

.995 

.960 

1000.0 

1.000 

.980 

5000.0 

1.000 

1.000 

TABLE  (4-2B) 

TROPOSPHERIC  ERROR  ADJUSTMENT 
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Elevation  (mr) 

er(Elev.)  (m) 

ee(Elev.)  (mr) 

0.000 

229.880 

.387 

17.453 

229.484 

.398 

34.907 

228.387 

.406 

69.813 

224.089 

.413 

104.720 

217.566 

.408 

139.626 

209.459 

.395 

174.533 

200.436 

.377 

261.799 

176.997 

.325 

349.066 

155.631 

.273 

523.599 

119.482 

.188 

698.132 

100.584 

.138 

872.665 

88.483 

.101 

1047.198 

79.248 

.071 

1221.730 

73.762 

.045 

1396.263 

71.110 

.021 

1483.530 

68.885 

.010 

TABLE  (4-3A) 


IONOSPHERIC  (DAY)  ERRORS 
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Altitude  (km) 

£x(Alt.) 

f2 (Alt.) 

100.0 

0.000 

0.000 

150.0 

.030 

.120 

200.0 

.090 

.230 

300.0 

.390 

.620 

400.0 

.750 

.950 

500.0 

.910 

1.000 

700.0 

.990 

.960 

1000.0 

1.000 

.890 

2000.0 

1.000 

.780 

5000.0 

1.000 

.600 

TABLE  (4-3B) 

IONOSPHERIC  (DAY)  ADJUSTMENT 
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Elevation  (mr) 

er(Elev.)  (m) 

ee(Elev.)  (mr) 

0.000 

64.587 

.116 

17.453 

64.465 

.120 

34.907 

64.130 

.122 

69.813 

62.850 

.126 

104.720 

60.899 

.125 

139.626 

58.491 

.122 

174.533 

55.778 

.117 

261.799 

48.859 

.101 

349.066 

42.642 

.085 

523.599 

34.138 

.062 

698.132 

27.432 

.045 

872.665 

23.805 

.030 

1047.198 

21.336 

.021 

1221.730 

20.117 

.013 

1396.263 

19.050 

.006 

1483.530 

18.288 

.003 

TABLE  (4-4A) 

IONOSPHERIC  (NIGHT)  ERRORS 
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Altitude  (km) 

f x (Alt . ) 

f 2 (Alt . ) 

100.0 

0.000 

0.000 

150.0 

.005 

.030 

200.0 

.090 

.170 

300.0 

.590 

.820 

400.0 

.860 

.990 

500.0 

.950 

.990 

700.0 

.995 

.920 

1000.0 

1.000 

.840 

2000.0 

1.000 

.760 

5000.0 

1.000 

.710 

TABLE  (4-4B) 

IONOSPHERIC  (NIGHT)  ADJUSTMENT 
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ATMOSPHERIC  £.  Imrod)  ATMOSPHERIC 


Fig.  4.1.  Tropospheric  induced  Range  Error,  er,  and  Elevation 
Error,  ee,  are  shown  as  a function  of  elevation  angle  in  (a)  and 

(c)  . The  fractional  adjustment  for  altitude  is  shown  in  (b)  and 

(d) . 
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Fig.  4.2.  Ionospheric  Daytime  induced  Range  Error,  er,  and 
Elevation  Error,  e^,  are  shown  as  a function  of  elevation  angle 
in  (a)  and  (c) . The  fractional  adjustment  for  altitude  is  shown 
in  (b)  and  (d) . 
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ATMOSPHERIC  Em  (mrocS)  ATMOSPHERIC 


Fig.  4.3.  Ionospheric  Nighttime  induced  Range  Error,  er,  and 
Elevation  Error,  ee,  are  shown  as  a function  of  elevation  angle 
in  (a)  and  (c) . The  fractional  adjustment  for  altitude  is  shown 
in  (b)  and  (d) . 
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4.3  Residual  Error 


In  general,  the  error  introduced  in  measured  target  position  as  a result 
of  refractive  bending  is  random  since  the  refractive  structure  of  the  atmosphere 
is  not  known  and  time  varying.  We  assume  that  some  statistical  method,  [9] , 

[10],  is  used  to  correct  for  the  effects  of  refraction  which  in  turn  produces  a 
residual  error.  The  resulting  residual  error  is  simulated  in  the  ARIES  program 
by  generating  a random  number  (for  each  coordinate  element  effected  by  refrac- 
tion) from  a distribution  whose  variance  is  related  to  the  uncertainty  in  the 
error  after  the  correction  method  is  applied.  The  random  numbers  for  each  co- 
ordinate element  are  used  in  the  measurement  equations  as  described  in  Section  2. 

The  residual  error  is  modeled  using  the  following  assumptions: 

1.  The  residual  error  is  Gaussian. 

2.  The  refraction  error  obtained  by  table  look  up  as  described  in 
Section  4.2  is  the  mean  refraction  error. 

3.  After  application  of  a refraction  correction  to  the  radar  position 
data,  the  standard  deviation  of  the  residual  error  is  assumed  to  be 
some  fraction  of  the  absolute  value  of  the  mean  refraction  error. 

4.  We  neglect  azimuthal  variations  on  the  refraction  terms.  Azimuthal 
refraction  arises  whenever  there  are  horizontal  gradients  of  refrac- 
tive index.  Studies  show,  [11],  that  such  horizontal  variations  are 
small  compared  to  the  vertical  variations  of  the  index  of  refraction. 
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4.3.1  Tropospheric  Residual  Error 

Based  on  the  assumptions  in  Section  4.3,  the  tropospheric  residual  error 
is  modeled  as  a zero  mean  Gaussian  process.  For  the  i— measurement , the 
standard  deviation  of  the  tropospheric  residual  errors  in  RAE  coordinates  for 
range  and  elevation  respectively  are 


°RTp^i)  P>p|  ARp(i)  | 

°ETp^  = pTlAEr(i) 


(4-10) 


where  Pp  is  the  fraction  representing  the  uncorrectable  part  of  the  error  (an 
ARIES  input  parameter)  and  ARp(i) , AEp(i)  are  obtained  from  Eqn.  (4-9)  after 
using  the  tropospheric  tables.  Azimuthal  refraction  errors  are  neglected.  In 
RUV  coordinates,  the  standard  deviation  of  the  residual  errors  for  the  i— 
measurement  are 


aj^Pp(i)  Prp|  ARp(i)  | 

Cjjpp(i)  - Pp | sin(E)  sin (A-AZ)  AEp ( i)  | 

a-vTp(i)  = Pjl  [cos  (E)  cos  (EL)  + sin(E)  sin(EL)cos(A-AZ)  ]AEp(i) 


(4-11) 


where 

E = target  elevation  angle  of  the  i—  measurement 

th 

A = target  azimuth  angle  of  the  i—  measurement 
EL  = elevation  angle  of  the  phased  array  broadside 
AZ  = azimuth  angle  of  the  phased  array  broadside 
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4.3.2  Ionospheric  Residual  Error 

Unlike  the  tropospheric  model,  the  ionospheric  model  is  a function  of 


radar  frequency,  as  can  be  seen  by  Eqn.  (4-6) . Tables  4-3  and  4-4  were 
computed  based  on  a radar  frequency  of  400  MHz.  For  a different  radar 
frequency,  the  errors  (ARj  and  AEj)  obtained  from  the  tables  using  Eqn.  (4-9) 
must  be  scaled  according  to  the  expression 


6Rj  = ARj (Fq/F) 2 

«Ej  = AEj (Fq/F) 2 


(4-12) 


where 

Fq  = frequency  used  in  generating  tables  4-3  and  4-4  (400  MHz), 
F = radar  frequency  used  in  simulation, 

I = subscript  to  indicate  ionospheric  errors. 


Based  on  the  assumptions  in  Section  4.3,  the  ionospheric  residual  error 
is  modeled  as  a zero  mean  Gaussian  process.  For  the  i— measurement,  the 
standard  deviation  of  the  ionospheric  residual  errors  in  RAE  coordinates  for 
range  and  elevation  respectively  are 


Opj ( i)  Pj | 6Rj ( i) | 

°Ei (i)  = Pj I (i) | 


(4-13) 


where  Pj  is  the  fraction  representing  the  uncorrec table  part  of  the  error  (an 
ARIES  input  parameter)  and  6Rj(i),  6Ej(i)  are  obtained  from  Eqns.  (4-9)  and 
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(4-12).  Azimuthal  refraction  errors  are  neglected.  In  RUV  coordinates  the 

th 

standard  deviation  of  the  residual  errors  for  the  i—  measurement  are 


°RI C i)  Pj|6Rj(i)| 

oui(i)  = Pj |sin(E)sin(A-AZ)6Ej(i) | 

Oyj(i)  = Pj | [cos (E) cos (EL)  + sin(E)sin(EL)cos(A-AZ)]6Ej(i) | 


(4-14) 


where  all  parameters  have  been  defined  in  Section  4.3.1. 
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5.0  IONOSPHERIC  SCINTILLATION 


The  ionospheric  scintillation  model  was  developed  by  R.  H.  Wand  [12]  from 
UHF  measurements  made  at  M.I.T.  Lincoln  Laboratory,  Millstone  Radar  Site,  over 
the  period  from  1971  to  1973.  Two  forms  of  the  model  are  considered,  a "0.1% 
probability  model"  and  a "worst  case  ever  model".  For  the  first  model,  the 
level  of  scintillation  specified  could  be  expected  to  be  exceeded  only  about 
.1%  of  the  time  for  a radar  operating  in  or  near  the  auroral  zone.  The  second 
model  is  based  on  measurements  made  at  Millstone  Hill  during  a period  of 
severe  magnetic  disturbances  (Kp  = 8+)  which  were  the  largest  fluctuations 
ever  observed  during  the  Millstone  satellite  tracking  program.  As  a comparison 
the  worst  case  ever  model  represents  a level  of  scintillation  approximately 
four  times  larger  than  the  .1%  probability  model. 

5 . 1 Ionospheric  Scintillation  Model 

Theoretical  and  experimental  considerations  indicate  that  the  amplitude 
fluctuations  produced  by  ionospheric  scintillation  may  be  represented  as  a 
function  of  samples  drawn  from  a Nakagami  M-distribution  [13]: 

7vM  C2M-1  -MB2  B > 0 

fM®  = ~m  - M > * C5-1) 

where  the  parameter  M = (2 and  Oq,  the  rms  amplitude  fluctuation,  is 
obtained  from  the  orbital  geometry  in  a manner  which  will  be  explained  shortly. 
Then  if  S’  represents  the  mean  of  the  signal  amplitude  (i.e.,  the  amplitude  if 
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no  scintillation  were  present) , and  if  A represents  the  received  signal 
amplitude  with  scintillation,  we  have 

A = G2  X (5-2) 

where 

G = (1  + aG)B2  (5-3) 

and  B is  drawn  from  the  distribution  in  (5-1) . The  G term  in  Eqn.  (5-2)  uhich 
represents  the  gain  or  attenuation  in  amplitude  due  to  the  scintillation  is 
squared  because  the  radar  measurement  is  a two-way  transmission. 

We  obtain  a^,  the  rms  amplitude  fluctuation,  as: 

°G  = sec  (5-4) 

where 


a 


0 


.04  for  the  "0.11  probability  model" 
.16  for  the  "worst  case  ever  model" 


and 


z 


zlz2 

zl+z2 


(5-5) 


where  Zp  z~,  and  h are  defined  in  Fig.  5-1.  The  zenith  angle  i at  the 
irregularities  is  obtained  from 
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1 


(5-6) 


. -1 
= sm 


cosE  \ 
^ ) 


where  rQ  is  the  earth's  radius.  Note  that  the  above  equation  applies  only  for 
z2  > °G  sh°u-*-^  t>e  made  zero  once  the  RV  target  penetrates  the  layer  of 

irregularities.  Also,  if  Eqn.  (5-5)  yields  a value  of  > .5,  then  we  should 
set  = .5. 

A given  value  of  B in  (5-2)  is  held  for  a correlation  time  x which  corre- 
sponds roughly  to  the  time  taken  for  the  penetration  point  at  height  h to 
traverse  a "typical"  ionospheric  irregularity.  In  general,  a mixture  of 
irregularity  scale  sizes  are  present,  but  the  size  which  is  most  prominent  in 
producing  the  amplitude  scintillations  has  a radius  of  one  Fresnel  zone,  \f\T, 
where  X is  the  observing  wavelength.  Then 


X 


(5-7) 


where  the  effective  penetration  point  speed  V is  given  by 


Vp  " IX(zl  * *2^  +f  >■ 


(5-8) 


• • 

and  E and  T are  the  rate  of  change  of  elevation  and  traverse  angle  respectively 
in  radians/sec.  If  the  RV  is  moving  in  such  a manner  that  is  small,  then  a 
minimum  value  of  100  m/sec  is  assumed,  corresponding  to  the  speed  of  the  most 
probable  size  irregularity. 
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Fig.  5.1.  Ionospheric  Scintillation  Model . 
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For  the  effects  of  the  ionospheric  scintillation  on  angle  measurements 
we  first  note  that  UHF  observations  at  Millstone  Hill  showed  that  angle 
scintillation  was  generally  much  more  pronounced  in  traverse  angle  than  in 
elevation  angle.  This  is  presumably  related  to  the  high  inclination  angle  of 
the  earth's  magnetic  field  and  the  fact  that  the  irregularities  are  elongated 
in  the  direction  of  the  field.  The  UHF  results  also  gave  approximately  a 
linear  relationship  between  rms  amplitude  scintillation  fluctuations  and  the 
rms  fluctuation  in  traverse  angle  Op. 

Op  = 40  Op,  mdeg  (5-9) 

Hence  Eqn.  (5-9)  provides  a means  of  predicting  a model  value  for  a,p  using  the 
model  value  of  a^.  The  fluctuations  in  elevation  angle  produced  by  ionospheric 
irregularities  will  be  neglected. 

The  probability  distribution  of  traverse  angle  fluctuations  may  be  taken 
as  Gaussian,  N(0,  a-p)  , while  the  correlation  time  may  be  assumed  to  be  the 
same  as  for  amplitude  fluctuations,  Eqn.  (5-7).  No  adjustment  to  the  traverse 
angle  fluctuation  is  necessary  to  allow  for  the  two-way  path  involved  in  a 
radar  situation. 

5 . 2 Generation  of  Nakagami-M  Random  Numbers 

The  Nakagami-M-distributed  random  variable  B of  Eqn.  (5-1)  is  generated  in 
subroutine  NAKGMI  by  one  of  three  possible  methods  depending  on  the  value  of  M. 
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For  M > 15  a Gaussian  approximation  to  the  Nakagami-M  distribution  is  used 
which  is  of  the  form 


B = {1  + 1/9M  + (l/SM)^}3/2 


where  RN  is  an  N(0,1)  Gaussian  random  variable. 

For  3 < M < 15  a transformation  of  uniform  random  variables  is  used.  Let 
{ly  be  independent  uniformly  distributed  random  variables  over  the  interval 
(0,1) . Then  for  M an  integer, 


B.  „ 
int 


(5-10) 


will  be  distributed  as  a Nakagami-M  random  variable.  If  M is  not  an  integer 

we  form  M.  . as 
int 


M 


int 


l Ml 


where  1*1  represents  the  closest  integer  value  and  then  form  B as 


B = 


Bint 


-a 


where  B^nt  is  generated  by  setting  M=M^nt  in  equation  (5-10)  and  where 


b = 


r CMint  + '±  / 

^int^^intV 


1 


r (M  + h) 

vC1  r(M) 
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a 


r(M. 


int 


+ h) 


i 


- b 


r (M) 
Hh  (M) 


The  technique  works  by  transforming  the  distribution  of  the  closest  integer, 

^int’  so  ^ias  same  mean  and  variance  as  the  desired  Nakagami-M 

2 

distribution.  This  method  successfully  passed  a x -test  on  its  density  func- 
tion for  values  of  M between  2 and  35.  As  a safety  precaution  this  technique 
is  not  used  for  M < 3 in  ARIES.  For  M > 15  the  Gaussian  approximation  is  used 
because  it  requires  less  execution  time. 

For  . 5 s M < 3a  probability  integral  transformation  [14]  is  used.  Let 

A X 

f(x)  = / fM(iOde  (5-n) 

o 

be  the  distribution  function  for  the  Nakagami-M  density  function  of  Eqn.  (5-1). 
Because  Y = F(X)  is  a monotone  increasing  function  for  X > 0 then 

X = F"1 (Y)  0 s Y $ 1 (5-12) 

where  F ^ is  used  to  denote  the  inverse  of  F.  If  we  select  a Y which  is 
uniformly  distributed  from  0 to  1 , then  by  the  probability  integral  trans- 
formation, X will  be  seen  to  be  a Nakagami-M  distributed  random  variable.  Al- 
though this  method  could  also  be  used  for  M > 3,  the  former  techniques  are 
computationally  much  more  efficient.  The  first  method  executes  about  20  times 
faster,  while  the  second  method  executed  from  7 to  20  times  faster  dependent 
on  the  value  for  M. 
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5.3  Effect  of  Ionospheric  Scintillation  on  Signal  to  Noise  Ratio 


The  value  for  B obtained  from  Subroutine  NAKGMI  is  used  in  Eqn.  (5-3)  to 
determine  G which  is  in  turn  used  to  modify  the  signal  to  noise  ratio  (SNR)  as 
shown  in  Eqns.  (3-9)  and  (3-10)  of  Section  3.2.1.  Because  SNR  is  a function 
of  power  rather  than  amplitude,  a term  appears  in  place  of  the  G^  term  of 
Eqn.  (5-2) . 

5.4  Calculation  of  Traverse  Angle  Fluctuations 

The  traverse  angle  fluctuations  are  obtained  by  combining  the  output  of 

subroutine  Gauss,  which  generates  an  N(0,1)  Gaussian  random  number,  and 

Eqn.  (5-9)  for  a^,.  We  then  obtain  the  traverse  angle  fluctuation  (in  radians), 

A™  , as 
Trav’ 


4Trav  * -* * * 04  °G  ^ 013> 

■ .6981  x 10'3  RN(t)  (5-14) 

where  RN(x)  is  an  N(0,1)  Gaussian  random  number,  and  the  t is  used  to  denote 
that  the  particular  random  number  sample  is  held  for  x seconds  before  being 
replaced  by  another  sample. 

The  azimuth  fluctuation,  AAg^nt,  becomes 

^Sint  = (5'15) 

where  E is  the  elevation  of  the  target.  If  the  radar  is  a phased  array,  then 
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the  azimuth  angle  fluctuations  are  transformed  into  fluctuations  in  U and  V 
coordinates  by  the  transformation 


4USint  ’ -cos®  coS(A-AZ)4A5int 
AVgint  = cos(E)  sin(A-AZ)  sin(EL) AAg^nt 


where  AUg^nt  and  AV 


EL  are  the  bores ite 


are  the  fluctuations  in  U and  V,  respectively,  AZ  and 
azimuth  and  elevation  directions  of  the  phased  array,  and 


A is  the  azimuth  of  the  target. 
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6.0  MULTIPATH  ERRORS 


Multipath  errors  are  considered  to  be  a result  of  two  interfering  paths 
for  the  return  from  a target.  The  first  path  is  a direct  return  to  the  main- 
lobe  of  the  beam.  The  second  path  is  an  indirect  return  due  to  reflection 
from  the  earth's  surface  into  a sidelobe  of  the  beam  (see  Fig.  6-1).  Note 
that  this  formulation  is  not  intended  to  model  multipath  errors  in  the  main 
beam  at  low  elevation  angles.  The  latter  type  of  multipath  gives  rise  to  much 
greater  error  magnitudes. 


6. 1 Multipath  Sidelobe  Level  Model 

The  sidelobe  level  multipath  error  is  primarily  an  angle  error.  The 
total  amplitude  of  the  return  is  insignificantly  affected  by  the  sidelobe 
portion  of  the  incoming  signal,  which  is  severely  attenuated  by  the  sidelobe 
level  and  ground  reflection.  However,  this  small  signal  is  still  sufficient 
to  cause  substantial  pointing  errors  in  the  angle  tracking  loop. 

A detailed  treatment  of  the  multipath  problem  is  available  in  Barton  [15] 
from  which  the  following  model  is  obtained.  We  start  by  assuming  a monopulse 
radar  pointing  at  a target  with  elevation  angle  G.  . If  the  sum  and  difference 
channels  are  denoted  by  E and  A,  respectively , then  the  effect  of  the  reflected 
signal  is  to  generate  a spurious  output  in  the  A-channel  which  gives  rise  to 
an  rms  elevation  error  of 


. p Qbw 

w k yrsu 

m 


(6-1) 
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TARGET 


18-3-1684  7 


RADAR 

ANTENNA 


EARTH  S SURFACE 


V^7^W7777/7777/77?Z;^?77?7Z 


Fig.  6.1.  Multipath  errors  due  to  direct  and  reflected  returns  from 
the  target. 
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where  0^w  is  the  3 db  elevation  mainlobe  beamwidth,  p the  reflection  coeffi- 
cient at  the  reflection  point  on  the  earth's  surface,  and  the  normalized 
monopulse  A pattern  slope.  SLL  is  the  Z mainlobe  to  A sidelobe  voltage  ratio 
averaged  over  the  first  Fresnel  zone  nominally  centered  20t  below  the  mainlobe 
axis  [see  ref.  [15],  figure  2]. 

The  program  obtains  SLL  from  a pre- stored  table  which  depends  on  the 
particular  type  of  radar  which  is  being  simulated,  p is  calculated  as 


P = | R±|p  D (6-2) 

where 

- is  the  reflection  coefficient  from  a flat  surface  for  the 
received  polarization  given  by  equations  (6-3)  to  (6-5)  . 

D - is  the  divergence  to  account  for  the  earth's  curvature. 

p - is  the  surface  roughness  factor.  It  is  given  by  p=exp{ -4irAsinK} 
where  A is  the  mean  surface  roughness  in  wavelengths. 

E - elevation  angle. 

The  reflection  coefficient  for  the  various  types  of  polarizations  are  as 
follows : 


vertical 


e sinE  - \ 

e 

27 

- cos  E 

e sinE  + \ 

e 

- cos^E 

(6-3) 


horizontal 


sinE  - \|  e - cos“E 


sinE  + e - cos^E 


(6-4) 


circular  |RC|  = h[\\\2  + |RJ|  - 2^  |Rh|cos($h-*v)  ] (6-5) 
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REFLECTION  COEFFICIENT,  |RC| 


Fig.  6.2.  Reflection  coefficient  for  circular  polarization,  1 R^. ! , as 
a function  of  elevation  angle  (Re(c)  = 10,  Im(e)  = .018). 
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where  e is  the  complex  permittivity  of  the  ground,  and  <j>  and  are  the  phase 

lags  on  the  reflection  for  vertical  and  horizontal  polarization,  respectively. 
Figure  6-2  shows  |R  | as  a function  of  elevation  angle  for  the  real  part  of 
the  permittivity  equal  to  10,  the  imaginary  part  equal  to  0.018,  and  a flat 
smooth  surface. 

If  we  denote  the  multipath  elevation  error  by  AE^^  then 

4EMul  ’ °EM  RN<4E)  (6‘6J 

where  RN  is  zero  mean  unit  variance  Gaussian  random  number  and  the  notation 
AE  indicates  that  a new  sample  is  obtained  when  the  elevation  angle  changes  by 
an  amount  AE.  The  value  for  AE  is  selected  as  half  the  mainlobe  width,  so 
that  these  errors  remain  correlated  as  long  as  the  multipath  error  is  due  to 
interference  from  a particular  sidelobe.  For  a new  sidelobe,  then  a new  ran- 
dom number  sample  is  selected. 

For  a phased  array  radar  the  elevation  angle  errors  must  be  transformed 
into  U and  V coordinates  by  the  following  expression 


AUMul  = sin(E)  sin(A-AZ)AENlul 


AV..,  = sin  (EL)  cos(A-AZ)  + cos(E)  sin  (EL)  ] AEj^ 


where  AU^-^  and  AV^^  are  multipath  errors  in  U and  V respectively,  AZ  and  EL 
are  the  bores ite  azimuth  and  elevation  directions  of  the  phased  array,  and  E 
the  elevation  of  the  target. 
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6.2  Calculation  of  Multipath  Errors 


The  reflection  coefficient  magnitude,  R,  and  the  rms  sidelobe  level  values 
SLL,  are  stored  in  ARIES  as  arrays  which  are  functions  of  the  elevation  angle. 
See  Tables  6-1  and  6-2. 

The  stored  sidelobe  level  table  is  based  on  a mainlobe  beamwidth  of  .5°. 

For  other  values  of  mainlobe  beamwidth,  BW,  we  define  a scaled  elevation  angle 
„ • 

E as 

f'  = E 

and  then  do  a SLL  table  lookup  using  E'  in  place  of  E.  For  any  value  of 
E*  > 90°,  the  sidelobe  level  table  entry  at  90°  is  used.  Note  also  that  the 
sidelobe  level  table  (and  this  model)  are  applicable  only  when  the  elevation 
angle  is  greater  than  twice  the  mainlobe  beamwidth.  The  reflection  coefficient 
used  is  identical  to  the  curve  in  Fig.  6-1  since  A=0,  D=1  are  assumed  for  the 
stored  tables. 

In  operation  table  look-ups  are  made  to  obtain  both  p and  the  rms  side- 
lobe level.  These  are  combined  with  and  0^  to  form  the  a™  of  Eqn.  (6-1). 
Then  the  value  of  obtained  is  multiplied  by  the  same  random  sample  used  on 
the  previous  multipath  calculation  unless  the  elevation  angle  has  changed  by 
half  a mainlobe  beamwidth,  in  which  case  a new  random  number  sample  is  selected 
The  value  selected  for  k depends  on  the  exact  type  of  radar  used.  For 
example,  a monopulse  radar  might  typically  have  a value  for  k^  = 1.44  and  thus 
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where 


EM 


R 0BW 
1.44  2 SLL 


.491  R 0 


BW 


SLL 


om  and  0g^  are  expressed  in  radians  by  ARIES. 
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Elevation  (mr) 

Reflection  Coefficient 

8.727 

.96883000 

34.907 

.88430000 

61.087 

.81100000 

87.266 

.74674000 

113.446 

.68984000 

139.626 

.63905000 

165.806 

.59337000 

191.986 

.55203000 

218.166 

.51440000 

244.346 

.47999000 

287.979 

.42867000 

331.613 

.38364000 

375.246 

.34380000 

418.879 

.30828000 

462.512 

.27642000 

506.145 

.24772000 

549.779 

.22176000 

593.412 

.19819000 

637.045 

.17675000 

680.678 

.15721000 

724.312 

.13937000 

794.125 

.11399000 

881.391 

.08702700 

968.658 

.06463000 

1055.924 

.04619800 

1134.464 

.03262000 

1221.730 

.02055900 

1308.997 

.01143200 

1396.263 

.00504200 

1483.530 

.00127450 

1570.796 

.00027296 

Table  6-1.  Reflection  coefficient,  R,  as  a function  of  elevation 
angle  (Re(e)  = 10,  Im(e)  = .018,  A = 0,  D = 1). 
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Elevation  (mr) 

Sidelobe  Level  (dB) 

17.453 

-22.42 

18.326 

-23.53 

19.199 

-24.61 

20.071 

-25.78 

20.944 

-26.92 

21.817 

-28.00 

34.907 

-28.00 

34.907 

-33.98 

104.720 

-33.98 

104.720 

-40.00 

174.533 

-40.00 

349.066 

-40.00 

523.599 

-40.00 

698.132 

-40.00 

872.665 

-40.00 

1047.198 

-40.00 

1570.796 

-40.00 

Table  6-2.  Tables  stored  by  ARIES  for  the  RMS  sidelobe 
level  table  as  a function  of  elevation  angle. 
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7.0  RADAR  CROSS  SECTION  COMPUTATION 


The  purpose  of  this  Section  is  to  present  a summary  of  the  radar  cross 
section  (RCS)  computation  performed  in  the  ARIES  Program. 

The  ARIES  program  gives  the  user  the  option  of  selecting  one  of  four 
methods  for  computing  the  RCS  of  each  target.  The  first  method  of  determining 
the  RCS  value  is  to  do  a table  look-up  on  a set  of  stored  static  pattern  tables 
in  which  each  table  represents  the  radar  static  pattern  of  a specified  object. 
The  other  three  methods  of  determining  the  RCS  are  based  on  mathematical 
models  which  represent  a cylinder  with  end  plates,  a sinusoidally  changing  RCS 
with  given  frequency  and  null  depth,  and  a sphere  (constant  RCS).  The  four 
options  provide  extreme  flexibility  in  the  type  of  analysis  which  can  be  per- 
formed with  the  ARIES  program. 

Two  of  the  above  options,  the  table  look-up  and  cylinder,  require  the 
radar  aspect  and  roll  angles  to  the  target  to  be  determined  before  the  RCS  can 
be  computed.  Since  the  target  motion  is  assumed  to  have  six  degrees  of 
freedom,  Euler's  equation  must  be  solved  to  determine  the  object's  rotational 
motion  about  its  center  of  mass.  A coordinate  transformation  is  then  performed 
which  permits  the  computation  of  radar  aspect  and  roll  angles  from  Eulerian 
angles,  target  state  vector  position  and  radar  site  location.  In  the  following 
sections  we  will  present  a summary  of  the  methods  used  to  compute  the  radar 
aspect  angle,  roll  angle  and  target  radar  cross  section. 
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Fig.  7.1.  Definition  of  Eulerian  Angles. 


Fig.  7.2.  Definition  of  Object  Body  Axis. 
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7.1  Euler ian  Angles 


The  Eulerian  angles  are  defined  by  Fig.  7-1  in  which  (x  y z)  is  obtained 
from  (X  Y Z)  by  a counterclockwise  rotation  first,  about  the  Z axis  by  an 
angle  <f>,  then  about  the  x axis  by  an  angle  0,  and  then  about  the  z axis  by  an 
angle  ip.  The  transformation  between  the  coordinate  set  X = (X,Y,Z)  and 
x = (x,y,z)  is  given  by 

x = AX  (7-1) 


and  conversely 

— 1 — T 

X = A x = A x 

where  the  matrix  A is  defined  by 


(7-2) 


A(o,<m) 


cos<pcos<p  -cosOs  in<(>  s inip 
-sirapcosp-cosQsinpcosip 
sinOsinp 


cosipsin<j>+cos0cos4>sinip  sinipsinO 

-sin<psin<p+cosOcos<pcosip  cosipsinO 
-sin0cosp  cosO  _ 


(7-3) 


In  the  application  considered  here,  x denotes  the  principal  body  axis  of  the 

-1  T 

rotating  body  and  X is  an  inertial  frame.  Note  that  A = A since  the  matrix 
A is  orthogonal. 


7 . 2 Coordinate  Systems 

The  program  uses  several  coordinate  systems  in  which  different  functions 
are  performed  or  measured.  In  this  section  we  will  establish  relationships 
among  all  the  coordinate  systems  used  in  the  RCS  portion  of  the  program. 
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7.2.1  Principal  Body  Axes 

The  principal  body  axes  are  defined  as  a right-handed  system  whose  origin 
is  the  center  of  mass  of  the  object.  This  set  of  axes  is  attached  to  the 
object  and  assumed  to  be  rotating  with  it.  Definitions  of  the  object  body 
axis  and  rotational  parameters  are  given  in  Fig.  7-2. 

7.2.2  Local  Horizontal 

One  set  of  the  input  parameters  specified  in  the  ARIES  program  are  the 
Euler  angles  which  give  the  initial  orientation  of  the  body  axis  with  respect 
to  the  local  horizontal.  In  this  case,  the  local  horizontal  is  defined  as  a 
right-handed  system  with  the  Z-axis  along  the  velocity  vector  and  Y-axis 
normal  to  the  trajectory  plane. 

This  coordinate  frame  is  specified  using  the  inertial  state  vector  of 
the  object  at  the  initial  time  t = t . The  axes  of  the  local  horizontal  are 
obtained  as  follows: 

The  initial  velocity  vector  and  displacement  vector  are  normalized: 


I R|  = (Xj  + Yj  + l\)h 


hRX 

X 1 
1 

hRY 

= 

V,/|R| 

1 

tr 

7* 

1 M 

VIRl 

(7-5) 


where  the  subscript  I indicates  a component  in  the  ECI  frame.  Note  that  h„ 
and  h^  are  in  the  trajectory  plane. 

We  now  define  the  local  horizontal  Y direction  as  normal  to  the  trajectory 
plane. 


Ny  = 

— z 3 

h^Y  - hRy  hzz 

% 

= 

h.Z  hRY  ' hRZ  hZX 

hyz 

hZX  hRY  ' hRX  hZY 

Normalizing  Eqn.  (7-6)  we  obtain: 
|H|  = (H^  + + H^)*5 


hYX 

%/lHl 

hYY 

= 

%/|H| 

*Vz 

HyZ/ | H | 
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The  final  step  is  to  define  the  local  horizontal  X axis  as 


hx  = hz  X(-hy) 


hxx 

— — 

Nr  hZZ  hZY  hYZ 

hXY 

II 

hYZ  hZX  ' hZZ  hYX 

hXZ 

hyX  hZy  - hzX  hyy 

(7-8) 


Note  that  h^  is  a unit  vector  since  h„  and  hy  are  orthogonal. 

The  transformation  matrix  which  takes  the  initial  local  horizontal  into 
inertial  space  is  obtained  by  using  Eqns.  (7-4),  (7-7)  and  (7-8).  The  matrix 
is  as  follows: 


hXX 

hXY 

hXZ 

hyy 

]\z 

hZX 

hZY 

\u7 

LL 

(7-9) 


7 . 3 Rotational  Motion  of  a Rigid  Body 

The  rotational  motion  of  a rigid  body  with  no  external  torques  applied  is 
given  by  the  solution  of  the  homogeneous  form  of  Euler's  equations.  This 
solution  leads  directly  to  the  determination  of  the  time  varying  Euler  angles 
between  body  axis  of  the  object  and  an  inertial  frame.  In  this  section  we 
present  a summary  of  the  method  used  in  the  ARIES  Program  to  solve  for  the 
rotational  motion. 
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7.3.1  Equations  of  Motion 

The  absolute  motion  of  a body  is  described  in  terms  of  the  orientation 


of  the  principal  body  axes  (x,y,z)  with  respect  to  an  inertial  coordinate 

frame  (X,Y,Z).  If  (u>  , u>  , w ) are  the  components  of  angular  velocity  of  the 

x y z 

body  in  the  (x,y,z)  frame,  then  the  angular  velocity  vector,  and  the 
momentum  vector,  M,  in  the  (x,y,z)  frame  can  be  written  as 


0) 

"i  b>  " 

X 

X X 

n = 

0) 

M = 

I U 

y 

y y 

0) 

Tw, 

Z 

z z 

• — 

where  1^,  I and  I_  are  the  principal  moments  of  inertia  of  the  body.  For  the 
exoatmospheric  portion  of  the  trajectory,  it  is  assumed  that  there  are  no 
external  torques  applied  to  the  body.  This  implies  that  angular  momentum  is 
conserved,  that  is 


dM\ 

-j=  1 + Q x M = 0 

dt  W) 


(7-10) 


which  is  the  vector  form  of  Euler's  equations.  In  component  form  Eqn.  (7-10) 
becomes 
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• - VM 

“x  - - hHJ  “y“z  X 


• -/VM 

y \ i . / x z 


y 


“z  - (W  “x“y 


/ 


(7-11) 


Two  integrals  of  Euler's  equations  are  known  from  the  laws  of  conservation  of 
energy  and  angular  momentum: 


_ 2 

T 2 

2 

Id)  + 

I to  + 

I d) 

= 2E 

X X 

y y 

z z 

t2  2 

t2  2 

t2  2 

.,2 

Id)  + 

I to  + 

I d) 

= M 

X X 

y y 

z z 

where 

2EI  < M2  < 2EI  for  I < I < I 
x z x y z 

or  where 

2EI  > M2  > 2EI  for  I > I > I 

X Z X y Z 

in  which  the  energy  E and  the  magnitude  M of  the  angular  momentum  are  given 
constants. 

2 

Manipulation  of  Euler's  equations  and  the  equations  for  E and  M will 

give  uncoupled  equations  for  (to  , to  , to  ) which  may  be  solved  in  terms  of 

x y z 

Jacobian  elliptic  functions.  In  some  cases  a more  trivial  solution  may  be 
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obtained  depending  upon  the  relative  magnitudes  of  the  principal  moments  of 

2 2 

inertia  (Ix,  I , I,)  and  whether  M > 2EI^  or  M < 2EIy.  For  the  details  of 
the  solution  the  reader  is  referred  to  references  [16]  and  [17]. 

7.3.2  Motion  in  Space 

Using  the  calculated  angular  velocity  (to  , u>  , to  ) obtained  from  the 
solution  of  Euler's  equation,  the  motion  of  the  body  with  respect  to  an 
inertial  frame  can  be  determined.  This  implies  that  we  must  determine  the 
time  varying  Euler  angles  between  the  body  axis  and  the  inertial  frame.  To 
simplify  the  resulting  equations  one  axis  of  the  inertial  system  is  selected 
to  be  along  the  invariant  angular  momentum  vector  M.  The  inertial  system  is 
therefore  defined  to  be  the  momentum  frame  taking  the  Z-axis  in  the  direction 
of  the  constant  M vector. 

The  components  of  angular  velocity  can  be  expressed  in  terms  of  the 
Eulerian  angles  and  their  derivatives,  as 

wx  = sin0M  sin*M  + °M  C0S*M  ) 
y = *M  sin0M  C0S*M  ’ °M  sin*M  (7_12) 

z = *M  COS0M  + *M  ) 

where  the  subscript  M denotes  that  the  angles  are  taken  with  respect  the 
inertial  momentum  frame. 

The  solution  of  Eulerian  angles  (0w,  $w,  follows  using  the  angular 
momentum  equations.  The  polar  angle  and  azimuth  of  the  momentum  Z-axis  with 
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respect  to  the  body  axes  x,  y,  z are  respectively  0^  and  Jsn-^  and  we  obtain 
on  taking  the  components  of  M along  the  axes  x,  y,  z the  expressions: 


|M|  sinG^  sin^M  = Ixu>x(t) 
|M|  sinOM  cos*M  = I wy(t) 
|M|  cosO^,  = iz«z(t) . 


(7-13) 


Hence 


cos0M(t)  = Izwz(t)/|M| 
tan*M(t)  = Vx(t)/yy(t) 


(7-14) 


which  give  the  angles  ^ and  as  a function  time. 

The  angle  <fy  does  not  appear  in  Eqn.  (7-13)  and  to  calculate  it  use 
Eqn.  (7-12),  which  expresses  the  components  of  w in  terms  of  the  time  deriva- 
tives of  the  Eulerian  angles.  Eliminating  from  Eqn.  (7-12),  then  using 
Eqn.  (7-13)  we  obtain 


|m|  (i  u)2  + i w2) 

1 1 v x x y y 

(I  U)  + I 0)  ) 
v x x y y 


(7-15) 


The  function  ^(t)  is  obtained  by  integrating  Eqn.  (7-15). 


7.3.3  Initial  Euler  Angles  to  Momentum  Frame 

From  the  foregoing  result,  define  the  initial  Euler  angles  between  the 
momentum  frame  and  the  body  axes  at  t as 


69 


(7-16) 


W = tan'1  [lx»x(t0)/I/y(to)] 

W = 0. 

The  transformation  from  body  axes  to  the  momentum  frame  at  t is  given  by 

x = AT[eM(t0),  *M(t0),  *M(t0)]x  (7-17) 

where  A is  the  matrix  defined  in  Eqn.  (7-3). 

7.4  Initial  Conditions 

The  following  initial  parameters  are  required  in  order  to  define  a solu- 
tion to  Euler's  equations.  The  parameters  are  the  three  rotational  rates 
(w  , a)  , id  ) , the  three  moments  of  inertia  (I  , I , I ) around  the  principal 

x y Z-  X y Z» 

body  axes,  and  three  Euler  angles  (e^,  tji^)  which  relate  the  initial  body 
axis  to  the  local  horizontal.  The  initial  Euler  angles  between  the  momentum 
frame  and  the  principal  body  axes  are  given  by  Eqn.  (7-16). 

The  solution  to  Euler's  equation  is  obtained  using  the  method  outlined 
in  reference  [16] . 

7 . 5 Radar  Aspect  and  Roll  Angle 

From  the  foregoing  results  we  are  now  in  a position  to  define  an  overall 
transformation  between  the  principal  body  axes  and  inertial  space  from  which 
the  radar  aspect  and  roll  angle  can  be  computed.  The  transformation  is  shown 
below.  Using  Eqns.  (7-1),  (7-3),  (7-9)  and  (7-17)  we  obtain: 
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(7-18) 


X = V I [0^  ,1^  ] A[ Oj^(tQ)  ,<J>^(to)  ,<J^(t0)  ] A , fyj(t)  >^(t)  ]x 

X = Tx 


This  transformation  can  be  diagramed  as  follows: 


We  can  now  determine  the  radar  aspect  and  roll  angles  in  the  following 
manner.  The  aspect  angle  0 , shown  in  Fig.  7-3,  is  defined  as  the  angle  be- 

3. 

tween  the  radar  line  of  sight  and  the  re-entry  vehicle  (RV)  roll  or  z-azis. 
That  is 


0 

a 


-1 

cos 


Bsm  1 
_!  Bsm  I - 


(7-19) 


where  R<^  is  the  vector  from  the  radar  site  to  the  RV. 

Bsm  = " XS’  ym  ' Ys»  zm  " zs^ 


(7-20) 


71 


and 


(-XM’  YM’ 
(Xs>  Ys, 

2 = (^y>  Zv > 


Z^)  is  RV  coordinates  in  inertial  space. 

Zg)  is  radar  coordinates  in  inertial  space. 

Z-,)  is  the  position  of  the  roll  or  z-axis  of  the  object  in 
inertial  space  which  is  given  by 


Z = 


(7-21) 


The  roll  angle,  0 , shown  in  Fig.  7-3,  is  defined  by 


0 = tan 

r 


-1 


Bsm* 


(7-22) 


where 


I ■ (Yx,  Yy,  Yz) 


is  the  position  of  the  pitch  axis  in  the  inertial 
system,  that  is 


Y = 


(7-23) 


X = (X^,  Xy,  X^)  is  the  position  of  the  yaw  axis  in  the  inertial 
system,  that  is 


X = 


(7-24) 
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j PROJECTION  OF  RgM 
/ INTO  x-y  PLANE 


Definition  of  Aspect  Angle,  0a,  and  Roll  Angle,  0r. 


Fig.  7.3. 


The  quadrant  in  which  0 lies  is  determined  from  the  signs  of  the 
numerator  and  denominator  of  Eqn.  (7-22) , according  to  the  following  convention. 
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7.6  Radar  Cross  Section 


The  radar  cross  section  can  be  computed  by  one  of  four  methods.  The 
methods  will  be  described  here.  The  value  for  RCS  computed  by  any  one  of  these 
methods  is  used  in  Eqns . (3-8)  and  (3-9)  to  determine  the  signal  to  noise 
ratio  of  the  target. 

7.6.1  Static  Pattern 

Static  pattern  table  look  up  is  used  when  RCS  input  option  1 or  2*  is 
selected.  The  static  pattern  table  consists  of  measured  narrow-band  radar 
cross  section  data  with  two  orthogonal  polarization  components.  The  tables  are 
constructed  by  processing  the  measured  data  to  eliminate  points  which  can  be 
fitted  by  linear  interpolation  between  adjacent  points.  Maximum  table  size  is 
25  roll  cuts  with  400  aspect  points  per  roll  cut.  To  reduce  core  storage 
requirements  for  the  tables,  the  400  aspect  data  points  are  each  truncated  to 
12  binary  bits  in  accuracy  and  packed  5 points  per  60-bit  word.  In  this  way 
up  to  4 tables  can  be  kept  core  resident  during  a simulation  run. 

The  program  does  a static  pattern  table  look  up  for  each  specified  table 
type  using  the  aspect  and  roll  angles  computed  in  Section  7.5.  Since  many 
targets  are  symmetric  in  seme  way,  the  aspect  and  roll  angles  are  converted 
to  an  appropriate  interval  for  the  specified  table.  All  roll  and  aspect 
angles  are  first  transformed  to  the  interval 

*Input  option  2 is  not  currently  implemented. 
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(7-25) 


-180°  s 0„  •£  180° 
a 

0 < 0r  < 180° 


and  then  transformed  to  the  interval  for  the  symmetry  of  the  target  according 
to  table  7-1. 


Symmetry 

Code 

Range  of 
0a 

Range  of 
°r 

1 

(-180°,  +180°) 

(0,  180°) 

2 

(0°,  180°) 

(0,  90°) 

3 

(-180°,  +180°) 

(0,  90°) 

4 

(0°,  180°) 

(0,  45°) 

5 

(0°,  180°) 

(0,  40°) 

6 

(-180°,  +180°) 

(0,  45°) 

7 

(-180°,  +180°) 

(0) 

8 

(0,  90°) 

(0) 

TABLE  7-1 
SYMMETRY  RANGES 
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7.6.2  Cylinder  With  End  Plates 


A target  with  RCS  input  option  3 selected  is  modelled  as  a cylinder  with 
end  plates.  This  option  provides  the  capability  of  simulating  the  dynamic 
motion  of  a tank  without  the  need  for  a table  look  up.  The  radar  cross  section 
is  computed  using  the  following  expression 


rA  sin0. 


RCS  = 


sinfk  i cosG  ) 
v aJ 


~i2 


COSO 


+ rr cot“0  J-.“(2kr  sinO  ) 

3 i.  3 


(7-26) 


where 


r = cylinder  radius 
i = cylinder  length 
0 = aspect  angle 

Jj  = first- order  Bessel  function 
A = radar  transmitting  wavelength 
k = 2tt/A 


Note  that  the  first  term  in  Eqn.  (7-26)  gives  the  scattering  from  the  cylind- 
rical surface  while  the  second  gives  the  scattering  from  an  end  plate. 

The  RCS  for  both  left  and  right  circular  polarizations  are  set  equal  to 
RCS  in  Eqn.  (7-26).  A plot  of  Eqn.  (7-26)  as  a function  of  aspect  angle  for 
a given  l,  r and  A is  shown  in  Fig.  7-4. 
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0 25  50  75 


ASPECT  (deg) 

Fig.  7.4.  Radar  Cross  Section  of  cylinder  with  end  plate  where 
r = 4.115  cm,  1 = 20.894  cm,  X = .8611  cm. 
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7.6.3  Sinusoidal  RCS 


The  radar  cross  section  for  RCS  input  option  4 is  a sinusoidal  function 
of  time.  This  can  be  used  to  simulate  the  rotational  motion  of  symmetric 
objects  without  requiring  static  pattern  tables.  The  functional  form  of  the 
RCS  for  both  right  circular  and  left  circular  polarization  is  given  by  the 
expression  below: 

RCS  = |[1  + (1-Rf)sin(ut)  + Rf]  (7-27) 

where 

A = peak  RCS  amplitude 

Rr  = null  amplitude/A 

u)  = RCS  scintillation  rate  in  rad/sec. 

Note  that  input  parameters  for  this  option  enter  the  sinusoid  peak  and 
sinusoid  null  in  dbsm.  The  values  for  A and  R^  in  Eqn.  (7-27)  are  computed 
from  these  values. 

7.6.4  Spherical  Target 

RCS  option  5 specifies  that  the  radar  cross  section  is  a constant.  The 
RCS  amplitude  of  both  right  circular  and  left  circular  polarization  is  set  to 
the  value  specified  by  the  input  parameter. 
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7 . 7 RCS  Data  Tape 


An  Edit  option  in  the  ARIES  Program  permits  the  user  to  write  a magnetic 
tape  containing  the  measured  radar  cross  section  history  for  each  target  during 
every  track  segment.  The  measured  radar  cross  section  of  the  target  is  com- 
puted using  the  true  RCS,  an  ionospheric  scintillation  term  and  measured  target 
range.  This  value  effectively  represents  the  RCS  a radar  system  would  compute 
from  received  signal.  The  object  of  writing  a tape  of  this  type  is  to  permit 
offline  processing  and  analysis  of  target  signature  data. 

The  expression  used  in  computing  the  measured  right  and  left  circular 
polarization  RCS  is 

4 

RCSi  G4  (7-28) 

where 

R.  = measured  range  of  target 
Rp  = true  range  of  target 

RCSp  = radar  cross  section  of  target  computed  according  to  Section 
7.5  where  i=l  represents  right  circular  polarization  and  i=2 
left  circular  polarization 

G = ionospheric  scintillation  term  given  by  Eqn.  (5-3) 

The  complete  form  of  the  RCS  data  tape  is  given  in  Section  4.11.3  of 
reference  [1]. 
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